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The centrifugal force reversal and X-ray bursts 
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Abstract. Heyl (2000) made an interesting suggestion that the observed shifts in QPO frequency in type I X-ray 
bursts could be influenced by the same geometrical effect of strong gravity as the one that causes centrifugal force 
reversal discovered by Abramowicz and Lasota (1974). However, his main result contains a sign error. Here we 
derive the correct formula and conclude that constraints on the M(R) relation for neutron stars deduced from 
the rotational-modulation model of QPO frequency shifts are of no practical interest because the correct formula 
implies a weak condition R t > 1.3-Rs, where Rs is the Schwarzschild radius. We also argue against the relevance 
of the rotational-modulation model to the observed frequency modulations. 
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1. Introduction 

It is rather well established that type I X-ray bursts 
are caused by thermonuclear explosions in the material 
slowly accreted onto the surface of a neutron star. During 
the burst, X-ray emission is often quasi-periodic (e.g., 
Strohmayer, Zhang, & Swank 1997a), possibly as a re- 
sult of rotational modulation of the inhomogeneities in 
the burning matter. During the outburst the burning layer 
expands by a factor of about 10 (Joss 1978; Paczyhski 
1983) and according to Strohmayer et al. (1997b) the sub- 
sequent contraction accounts for the observed frequency 
shifts. Heyl (2000) noticed that in calculating these fre- 
quency shifts one should take into account various rela- 
tivistic effects such as the Lense-Thirring frame-dragging 
and the Abramowicz-Lasota (Abramowicz & Lasota 1974, 
1986; Barrabes et al. 1995 ) centrifugal-force reversal. 
Unfortunately Heyl used an unnecessarily complicated 
derivation of these effects and obtained an incorrect re- 
sult. Below we obtain the correct formula and discuss un- 
der what assumptions it can be applied, if at all, to QPOs 
frequency shifts observed during type I X-ray outbursts. 



In this case, one has a well defined problem to solve: to 
calculate the change in the angular velocity Ail, knowing 
the change in radius Ar (that corresponds to the thickness 
of the atmosphere) , and assuming that the specific angular 
momentum I is conserved: 



/ /=const 
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Because the spin of a neutron star in the rotational- 
modulation model may be assumed to be small, all cal- 
culations can be made in a metric in which only first or- 
der rotational corrections to the Schwarzschild metric of 
a non-rotating, spherical star are considered. With this 
accuracy, the relevant metric components are (Hartle and 
Thorne, 1968): 
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2. Frequency shifts 

Let us assume, as Heyl did, that as the burning region 
of the atmosphere (or a hot spot) expands and then con- 
tracts, the specific angular momentum of a fluid element 
that emits observed X-rays is conserved: I = —u v /ut, 
where u v , u t are components of the four- velocity. 
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Here M = GM*/c 2 is the mass of the star in geometrical 
units, and J = GJ*/c 3 is its angular momentum, also 
expressed in geometrical units. From now on, all formulae 
in this paper are given with the same accuracy to linear 
terms in rotation. It will be convenient to define the radius 
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of gyration f and the angular velocity of frame dragging 
u by 
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see e.g. Abramowicz, Miller and Stuchlfk (1993) for covari- 
ant definition and discussion of the meaning of the radius 
of gyration. 

Heyl (2000) attempted to solve the problem by using a 
method that was correct in principle, but by far unneces- 
sary complicated: he transformed the metric (2)-(4) to a 
reference frame that corotates with the star, derived and 
integrated in this frame the geodesic equation, discussed 
its solution introducing the Coriolis force, and transformed 
the solution back to the non-rotating frame. 

One may derive the desired result in only one line of 
calculations, realizing that since f2 = u v /u l , one gets from 
our definition of I and from Eq. (6) : 

n = ^ + u>, (7) 

to first order in rotational effects. 

^From this, just by direct differentiation (with I = 
const, 6 = const) one immediately arrives at, 

* a 0(,- 25£)-' . (8) 

ar r \ r J \ r ilr a J 

The sign of the (last) term proportional to J directly fol- 
lows from the monotonic decrease of uj with radius. 

It is remarkable, that to first order in stellar rotation, 
an initially uniformly rotating shell continues to rotate 
uniformly as it expands or contracts. Our formula shows 
that near r = 3M the shift in frequency Ail should be 
smaller than that predicted by Newtonian theory. This has 
correctly been noted by Heyl (2000), and it reflects the 
fact, known previously (Abramowicz & Prasanna 1990), 
that in the gravitational field of a non-rotating body, mat- 
ter which moves on nearly circular orbits with constant 
angular momentum experiences no shear (in the sense that 
dVt/dr — 0) exactly at the location of the circular photon 
orbit (r = 3M), i.e., at the location where the centrifugal 
force reverses. Eq. (^) shows that for matter moving on 
nearly circular orbits with constant angular momentum 
in the gravitational field of a (slowly) rotating body, zero 
shear occurs at some radius ro < 3M, contrary to Heyl's 
result, but in accordance, e.g., with the well known fact 
that the greater the angular momentum of the black hole, 
the smaller the radius of the (corotating) circular pho- 
ton orbit r p h- (Note, however, that only for non-rotating 
bodies ro = r p h, in general this is not true, because for 
rotating bodies r p h depends only on M and J, while r$ 
depends, in addition, on fL) 

In any case, the derivative of f2 in eq. (8) is positive 
only for r < 2.6M + [5//(2i? 2 ) — 1]. For realistic neutron 



star models, the moment of inertia is lower than that of 
a uniform Newtonian sphere, / < 2Mi? 2 /5, so the zero 
shear radius is less than the causality limit for neutron 
star radii, i?* > 2.8M (Haensel, Lasota & Zdunik 1999). 

3. Discussion 

We may ask if the formula (||) is relevant to the prob- 
lem of frequency shifts observed in QPOs during type I 
X-ray bursts and to the constraints on the mass-radius 
relation for neutron stars. For this, the Strohmayer et al. 
(1997b) model relating the QPO frequency to the neu- 
tron star's spin and the frequency shift to movements of 
the atmosphere must be correct, and during the expan- 
sion and contraction the specific angular momentum I of 
the QPO source should be conserved. We are skeptical on 
both counts. 

Muno et al. (2001), find that the presence of a ~ 600 
Hz quasiperiodic oscillation in the X-ray emission is corre- 
lated with radius expansion in X-ray bursts, but that the 
presence of a ~ 300 Hz QPO is not correlated. This in itself 
raises doubt as to the general validity of the rotational- 
modulation model of these QPOs. Further, we note that 
an independent confirmation of the rotational period of 
the neutron stars is needed before the frequency observed 
can be identified with the spin rate. Second, no long- 
lasting inhomogeneities have as yet been demonstrated in 
computations of the burst evolution. Third, even if the 
model is correct, there seem to be so many uncertainties 
involved in the motion of a hot spot in a stellar atmo- 
sphere, that we do not see how any frequency change can 
be reliably identified with the properties of the space-time 
metric. Suppose that no frequency change is observed. 
Would we claim that this is because the radius of the star, 
r = R coincides with the zero-shear radius ro, or should 
we rather suppose that the hot spot is stationary at the 
surface of the rotating star? Is the frequency decrease ob- 
served in 4U 1636-53 (Strohmayer 1999), due to R < r 
or is some other effect involved (perhaps the hot spot is 
rising like a balloon)? Finally, we know of no atmospheric 
phenomenon on Earth which would resemble the proposed 
model — castellanus clouds do not show a systematic west- 
erly bend with altitude, and large scale motions which 
can be ascribed to Coriolis forces, such as the jet stream 
or trade winds, are related to latitudinal motions of air 
masses, not vertical. 

Is I = —u^/ut conserved when the QPO is observed? In 
the case of axially symmetric, stationary motion of a per- 
fect fluid, both —u t (P + p)/n and u v (P + p)/n, where P, 
p and n are respectively the pressure, energy and particle 
densities, are conserved (see, e.g., Bardeen 1973) and be- 
cause I is simply their ratio, it is conserved as well. It is not 
clear what relevance this has to the frequency observed in 
X-ray bursts. After all, if the frequency is related to the 
stellar rotation it must be related to non-axisymmetry, 
breaking the basic assumption under which I is conserved 
for perfect fluids. 
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A naive reading of the suggestion that the X-ray burst 
frequency may be caused by a dark or bright fluid element 
conserving angular momentum as it changes altitude, sug- 
gests the motion of a hovercraft. This would be a particle 
free of azimuthal forces but supported vertically, and the 
relevant question would be the change of its azimuthal 
velocity or simply angular frequency, as seen by a distant 
observer, when the particle changes its radial position. The 
energy per unit rest mass — u t of the particle would not be 
conserved, but its angular momentum per unit rest mass 
u v = 1(1 — 2M/r) 1 / 2 would, and hence, in this case, 

dn _ CI ( 2M\ _1 
dr r \ r J 

could be a better model for frequency shift than (|J). 

In conclusion, it seems unlikely that Heyl's idea can 
ever be used to constrain the mass-radius relationship of 
a neutron star. Unfortunately, it seems that the effect two 
of us discovered 27 years ago has yet to be observed. 
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